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ABSTRACT 


In  this  thesis,  the  general  solution  of  the  biharmonic  equation 
in  Cartesian  coordinates  is  used  to  solve  the  following  problems: 

1.  Uniformly  loaded  rectangular  plates  with  all  edges  clamped. 

2.  Uniformly  loaded  rhombus  plates  with  all  edges  clamped. 

3.  Uniformly  and  hydrostatically  loaded  circular  segmental 
plates  with  all  edges  clamped. 

4.  Uniformly  loaded  circular  plates  with  diametrically  opposite 
flat  sides  and  all  edges  clamped. 

5.  Creeping  flow  in  rectangular  cavities. 

6.  Creeping  flow  in  triangular  cavities. 

7.  Creeping  flow  in  circular  segmental  cavities. 

The  general  solution  of  the  bi harmonic  equation  is  written  in 
terms  of  infinite  series  of  harmonic  and  bi harmonic  functions  which  are 
derived  by  using  the  complex-variable  method.  The  coefficients  of  the 
infinite  series  are  determined  by  point-matching  technique. 

The  numerical  results  have  been  compared  with  the  available  solu¬ 
tions  in  the  literature.  It  is  seen  that  the  numerical  results  of  the 
present  analysis  are  very  satisfactory. 

The  purpose  of  this  thesis  is  to  show  that  a  class  of  boundary- 
value  problems  can  be  approached  by  using  bi harmonic  functions  in 


Cartesian  coordinates. 
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NOTATIONS 

A 

Aspect  ratio  of  rectangular  cavity  (=H/a). 

A  ,B  ,C  ,D 
m  m  m  m 

Coefficients  of  the  general  solution  of  biharmonic 

equation  in  polar  coordinates. 

amn 

Coefficient  of  the  series  solution  of  the  bi harmonic 

equation  in  doubly  infinite  series. 

a 

Length  dimension  in  x-di recti  on. 

b 

Length  dimension  in  y-direction. 

bi 

Coefficient  of  the  bi harmonic  function  G . . 

ci 

Coefficient  of  the  harmonic  function  g . . 

D 

=  Et^/12(l-y^) ,  Flexural  rigidity  of  a  plate. 

E 

Modulus  of  Elasticity. 

Gi 

Bi harmonic  function. 

g . 

Harmonic  function. 

H 

Depth  of  the  rectangular  cavity. 

k,  £ 

Real  numbers. 

Mx*VMn 

Bending  moments  per  unit  length  of  sections  of  a  plate 

perpendicular  to  X  axis,  Y  axis  and  unit  normal, 

respectively. 

n 

Outward  normal  of  the  boundary. 

q.  q0 

Load  intensity. 

|3n,c'n’rn,sn 

Real  constants. 

M,  N,  P 

Positive  integers. 

R 

Reynolds  number. 

s 

Length  dimension  in  y-direction  or  along  the  boundary. 

Velocity  of  the  moving  wall  for  creeping  flow  problem. 
Velocity  of  the  creeping  flow  problem. 

Deflection  of  the  plate. 

Stream  function  of  the  creeping  flow  problem. 

Complex  constant. 

Complex  constant. 

Poisson's  ratio  (=0.3  for  all  the  numerical  results 
of  the  plate  problem). 


CHAPTER  I 


INTRODUCTION 


1.1  Plate  Equation  and  Its  Solutions 


The  governing  differential  equation  of  a  laterally  loaded  thin 
plate  subjected  to  small  deflections  is 


V4w 


94w 

9x^ 


+  2 


94w 


2  2 
9x  ay^ 


94w 

“T 

9y 


D 


In  the  case  of  a  clamped  edge  the  boundary  conditions  are 

r\.  . 

w  =  0  and  Jn  =  where  n  outward  normal. 


The  general  solution  of  differential  equation  (1)  can  be  written  as, 

w  =  w,  +  w^  , 
h  p 

where  wp  is  a  particular  integral  of  the  differential  equation  (1)  and 
w^  is  the  solution  of  the  homogeneous  bi harmonic  differential  equation 

4 

V  w,  =  0  . 

h 

The  general  series  solution  of  the  homogeneous  biharmonic  equation  in 
polar  coordinates  is  [1], 

wh  =  A0  +  V2  +  c0  lo9  r  +  DQr2  log  r 

+  ( A-j r  +  B^r3  +  C-j r“ 3  +  D-j r  log  r)  cos  0 


1 


2 


+  (A-jr  +  B-jr2  +  C-jr"^  +  D-jr  log  r)  sin  6 

+  1  ( Am^  +  B  r~m  +  Crm+2  +  D  r~m+2)  cos  me 

„_o  m  m  m  m  ' 

m=2 

+  I  (Amrm  +  Bmr’m  +  Cmrm+2  +  D'r“m+2)  S1'n  m0  (2) 

m  m  m  m 

m=2 

This  series  solution  (2)  has  been  successfully  used  by  Conway  [2], 

Cheng  [3],  Hulbert  [4],  Lo  [5]  and  many  others  to  solve  the  elasticity, 

plate  and  heat  transfer  problems  where  the  coefficients  A  .  B  .  C  , 

m  m  m 

Dm>  A^,  B^,  C^  and  were  determined  by  the  point  matching  technique. 

It  appears  that  the  general  solution  of  the  biharmonic  equation 
in  Cartesian  coordinates  was  first  given  by  Thorne  [6].  Sparrow  and 
Haji-Sheikh  [7]  set  up  the  general  solution  of  the  biharmonic  equation 
in  terms  of  G^1  and  g. ,  where  g.'s  and  G. 's  are  harmonic  and  biharmonic 
functions  respectively.  Sparrow  and  Haji-Sheikh  have  used  these 
functions  g.  and  G^  to  solve  the  flow  and  heat  transfer  in  ducts  with 
arbitrary  thermal  boundary  conditions.  In  this  thesis,  the  general 
solution  of  the  bi harmonic  equation  in  Cartesian  coordinates  will  be 
derived  by  complex  variable  method  and  polynomial  method.  It  is  noted 
that  the  functions  G. 's  obtained  are  not  completely  identical  with 
those  given  in  Reference  [7]. 


Derivation  of  the  biharmonic  function  G.  was  not  discussed  in  [7]. 
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1.2  Purpose  and  Problems 

The  general  solution  of  the  bi harmonic  equation  in  Cartesian  co¬ 
ordinates  was  first  used  successfully  by  Sparrow  and  Haji -Sheikh  [7] 
in  solving  laminar  forced  convection  heat  transfer  problem.  The  pur¬ 
pose  of  this  thesis  is  to  show  that  the  general  solution  of  bi harmonic 
equation  in  Cartesian  coordinates  can  be  used  to  solve  the  plate  and 
creeping  flow  problems  governed  by  bi harmonic  equations.  It  is  noted 
that  these  functions  G.  and  g^.  have  not  been  used  in  solving  these 
problems  in  the  past.  Approximate  solutions  are  found  by  employing 
point-matching  technique  to  obtain  the  coefficients  of  the  functions 
g.  and  G.  for  the  following  plate  and  creeping  flow  problems: 

1.  Uniformly  loaded  rectangular  plates  with  all  edges  clamped. 

2.  Uniformly  loaded  rhombus  plates  with  all  edges  clamped. 

3.  Uniformly  and  hydrostati cal ly  loaded  circular  segmental 
plates  with  all  edges  clamped. 

4.  Uniformly  loaded  circular  plates  with  diametrically  opposite 
flat  sides  and  all  edges  clamped. 

5.  Creeping  flow  in  rectangular  cavities. 

6.  Creeping  flow  in  triangular  cavities. 

7.  Creeping  flow  in  circular  segmental  cavities. 

It  is  noted  that  some  of  the  problems  mentioned  above  cannot  readily 
be  handled  by  using  the  general  solution  of  the  biharmonic  equation 
in  polar  coordinates.  Those  problems  are  clamped  uniformly  loaded 
rectangular  plate  with  large  aspect  ratio  and  clamped  uniformly  loaded 
rhombus  plate  with  small  angle  a. 


. 
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CHAPTER  II 


THEORETICAL  ANALYSIS 


2.1  Derivation  of  the  Bi harmonic  Function  g1  and  Harmonic  Function 
gl  by  Complex  Variable  Method 


The  differential  equation  (1)  in  complex  coordinate  system 
(z,  z)  can  be  written  as; 


where 


3z  3z^ 

z  =  x  +  iy,  7 
w(z,  I)  =  w(-2-g- 2 
q(z,  z)  =  q (-■ 


1 


3  (z,  z) 


1  =  /TT 
w(x,  y)  and 
q(x,  y)  . 


It  is  noted  that  w(z,  z)  and  q(z,  z)  are  both  real  functions.  The 

general  solution  of  equation  (3)  is  [8], 

» 

w(z,  7)  =  ffff  q(z,  7)  dz  dz  d7  d 7 

+  7  $(z)  +  z  7(7)  +  'F(z)  +  7(7)  . 


Here,  [z  $(z)  +  z  $(z)]  and  ['F(z)  +  'F(z)] 
are  bi harmonic  and  harmonic  functions,  respectively. 


1.  See  Appendix  1  for  forms  of  g.  and  G. . 


(3) 


4 


5 


Let  w(z,  z) 


where 


w,  +  w  , 
h  p 

■ygjj  ffff  q(z,  7)  dz  dz  dz”  dz" 
z  $(z)  +  z  J(J)  +  'i'(z)  +  ¥(7). 


( 3A) 


Here,  $(z)  and  'F(z)  are  analytic  functions  which  can  be  written  in 
the  following  forms; 


$(z)  = 

n 


Hz)  = 

n 


where  <J>  and  ip  are  unknown  complex  constants.  Now,  d>  and  ip  can  be 
n  n  r  n  n 

wri tten  as ; 

On  ’  Pn  + 

^  ■  rn  +  isn 

where  pn,  qn,  rn  and  sn  are  unknown  real  constants. 

Equation  (3A)  can  be  expanded  as. 


».(z,  z)  =  (x  -  iy)  l  (p  +  iqn)(x  +  iy)n 

n=l 


+  (x  +  iy)  l  (p  -  iq  )(x  -  iy)n 
n=l 

00  00 

+  l  (rn  +  isn)(x  +  iy)n  +  l  (rn  -  isn)(x  -  iy)n 

n=0  n=0 
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■  2Pn(x2  +  y2)  +  2p2(x3  +  y2x)  -  2q2(x2y  +  y3)  +  2p3(x4  -  y4) 

O  o  00 

-  4q3(x  y  +  xy  )  +  (x  -  iy)  £  (p  +  iq  ) (x  +  iy)n 

n=4  " 

00 

+  (X  +  iy)  l  (p  -  iq  )(x  -  iy)n  +  2r  +  2r  x  -  2s,y 

P=4  ii  II  ooi 

p  p  00 

+  2r2(x  -  y  )  -  4s2xy  +  £,  (r  +  is  )(x  +  iy)n 

L  n=3  n  n 

+J3  (rn  '  isn)(x  -  • 


Finally,  w^(z,  z)  can  be  written  as; 


wh(z,  z)  =  2p1G1  +  2p2G2  -  2q2G3  +  2p3G4  -  4q3G5 


’09l  +  2rl®2  *  2sl93  +  Zr2*4 


’n 


=  y  bG+y  eg 

n=l  "  n  n=l  n  n 


where  c  =  ir.  or  is. 

n  i  i 


bn  -  kp.  or  kq. 

and  k,  i  are  certain  real  numbers. 

The  first  five  terms  of  b^  and  c..  are 

b1  «  2p1 ,  b,  =  2p2,  b3  =  -q2,  b4  =  2p3,  bg  =  -4q3>  and 

c-j  ~  2r0,  c2  =  2 r* i  #  c3  ~  ~ 2s -j  *  C4  =  2r2»  -  ~2s2  . 
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2.2  Derivation  of  the  Functions  g.  and  by  Polynomial  Method 


The  functions  g.  and  can  also  be  obtained  by  using  a  doubly 
infinite  power  series. 

Suppose  is  regular  in  the  origin,  then  wh  can  be  written  as 
the  following  doubly  infinite  series. 


oo 


00 


l  l 

m=0  n=0 


mn 


xm  yn 


It  is  noted  that  the  coefficients  amn  are  not  all  linearly  independent. 

The  relationships  among  the  coefficients  a  „  can  be  found  after  satis- 

r  3  mn 

4 

fying  the  equation  V  wh  =  0. 

Let  w,  =  \  ?  a  xm  yn,  where  m  +  n  s  p  and  p  £  4. 

n  m=0  n=0  mn 

There  are  [(p  +  l)(p  +  2)/2]  coefficients  of  amn.  Further  the  bi¬ 
harmonic  equation  V4  w^  =  0  gives  rise  to  [(p-3)(p-2)/2]  "constraint" 
equations  which  define  the  relationships  among  the  coefficients  amn. 

The  total  number  of  independent  coefficients  amp  is 

(P  +  !)(P  1  1).  _  (P  ~  3^(P  ~J1  =  4p  -  2  (4A) 


The  equation  (4A)  is  valid  for  p  *  4  only.  Finally,  the  general  solu¬ 
tion  of  the  bi harmonic  equation  can  be  written  as; 


w 


h 


n  n  m  n 

-  I  I  amn  x  y 
m=0  n=0 


(2p-3)  (2p+l) 

l  b.  G.  +  l  c.  g. 

.1=1  J  J  1=1  1 


(4B) 
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where  m  +  n  <  p  and  p  <  4  and  the  coefficients  and  c..  are  now  all 
linearly  independent. 

As  p  -►  °°,  equation  (4B)  can  be  written  as; 


00  00 


w. 


1.  1.  amnx"V 


00 


00 


m=0  n=0 


mn 


1  k-  G.  +  I  ci  9-; 
j=l  J  J  i =1  1  1 


(4C) 


For  more  details,  see  Appendix  2. 


It  is  noted  that  the  sequence  xmyn  , 

fl  U  j  I  5  L  j  J  j  ,  #  , 


is  complete  approximating  basis  for  a  continuous  function  w(x,  y)  which 
is  finite  at  the  origin. 

The  functions  and  g^  may  be  regarded  as  complete  approximating 
basis  for  the  general  solution  of  the  biharmonic  equation  provided  the 
solution  is  regular  in  the  origin. 

A  computer  programme  was  written  to  obtain  results  using  the 
point-matching  technique.  The  Jordan  reduction  method  was  used  to 
solve  the  simultaneous  linear  algebraic  equations  to  obtain  the  coef¬ 
ficients  b.j  and  c.  of  the  general  solution  and  g^ . 


■ 
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CHAPTER  III 


APPLICATIONS  TO  SMALL  DEFLECTIONS  OF  LATERALLY 
LOADED  PLATES  WITH  ALL  EDGES  CLAMPED 

3.1  Uniformly  Loaded  Rectangular  Plates  with  All  Edges  Clamped 

The  problem  of  bending  of  clamped  rectangular  plates  under  uniform 
loading  has  attracted  many  investigators  in  the  past  and  various  approxi¬ 
mate  methods  of  solution  are  available  in  the  literature.  The  most 
notable  is  the  well-known  solution  of  Timoshenko  [1]  where  the  Levy 
solution  for  a  uniformly  loaded  simply  supported  rectangular  plate 
was  employed  and  edge  moments  were  applied  in  such  a  way  that  the  re¬ 
sulting  slopes  on  all  four  edges  are  approximately  zero.  This  method 

i 

of  solution  leads  to  a  doubly  infinite  set  of  simultaneous  equations. 

The  improved  method  of  solution  by  point-matching  was  discussed  by 
Conway  [2]  and  more  recently  by  Lo  and  Niedenfuhr  [9],  Due  to  the 
symmetry  of  the  problem  (see  Fig.  1(a)),  the  general  solution  can  be 
wri tten  as , 

w  =  D  +  x2y2  +  +  C191  +  blGl  +RI/C4k94k  +  b4kG4k^  ^ 

The  numerical  results  after  satisfying  the  boundary  conditions  w  =  0, 

~  =  0  at  5  equally  spaced  points  along  the  edge  AB  and  w  =  0,  ~  =  0 
at  4  equally  spaced  points  along  the  edge  BC  are  shown  in  Table  1  and 
compared  with  the  results  from  Timoshenko  and  Woinowsky-Krieger  [1]. 

It  is  seen  that  the  results  from  the  present  analysis  are  very 
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Table  1  Deflections  and  Bending  Moments  for  a  Uniformly  Loaded 
Rectangular  Plate  With  All  Edges  Clamped,  Poissons  Ratio  =0.3 
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satisfactory ,  Incidentally,  independent  calculation  shows  that  the 


wi th  —  =  1.0  in 
a 


value  0.0231  given  for  (M 


Timoshenko  and  Woinowsky-Krieger 's  book  [1]  is  slightly  in  error  and 
the  correct  value  is  0.0229. 

The  boundary  errors  for  deflection  and  slope  are  less  than  1.6% 
and  3.8%,  respecti vely ,  of  the  maximum  values  of  the  deflection  and 


slope  inside  the  plate  when  <  3.  In  general,  the  boundary  errors 
increase  as  the  ratio  —  increases.  It  is  noted  that  forms  of  particular 
solutions  such  as  w^  =  (x^  +  y^)  and  wp  =  ^  (x^y^)  may  be  used 


and  they  all  lead  to  the  same  results.  In  actual  computation,  it  is 
more  convenient  to  use 


3.2  Uniformly  Loaded  Clamped  Rhombus  Plates 

The  symmetry  of  the  problem  (see  Fig.  1(b))  shows  that  the  equa¬ 
tion  (5)  can  also  be  used  for  the  rhombus  plate  with  all  edges  clamped. 

The  numerical  results  after  satisfying  the  boundary  conditions  w  =  0, 

W 

=  0  at  9  equally  spaced  points  along  the  edge  AB  are  listed  in 
Table  2  together  with  the  available  results  from  the  literature  [1,10,11]. 
The  agreement  among  the  various  results  is  good  except  in  the  case 
a  =  10°  where  a  maximum  error  of  about  3%  was  observed  for  the  maximum 
deflection. 

The  edge  bending  moment  distribution  are  shown  in  Fig.  2.  For  the 
angle  a  <  20°,  some  irregularities  exist  around  the  acute  corner.  This 
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Table  2  Deflections  and  Bending  Moments  for  a 
Uniformly  Loaded  Rhombus  Plate  With 
All  Edges  Clamped,  Poisson's  Ratio  =  0.3 


Values  at 

center  of  plate 

Max.  value 

along  edge 

a 

p  n 

p 

p 

p 

At  distance 

degrees 

wxlO  /(qa h/D) 

MxlO/qa^ 

M  xlO/qa*1 

M  /qa 

may'  ^ 

s/( a/sin  a) 

y 

* 

1 II  CL  A 

from  obtuse  corner 

10 

3.021 

1.359 

0.6921 

-0.2677 

0.14 

0)* 

2.819 

1.347 

0.7200 

-0.2787 

— 

15 

2.436 

1.198 

0.6906 

-0.2499 

0.16 

20 

1.964 

1.054 

0.6766 

-0.2175 

0.18 

(1) 

1.961 

1.055 

0.6746 

-0.2095 

— 

(2)* 

1.974 

1 .056 

0.6771 

-0.2120 

0.18 

25 

1.567 

0.9190 

0.6527 

-0.1847 

0.22 

(2) 

1.568 

0.9188 

0.6528 

-0.1859 

0,22 

30 

1.229 

0.7914 

0.6172 

-0.1620 

0.31 

(1) 

1.230 

0.7916 

0.6176 

-0.1624 

— 

(2) 

1.230 

0.7916 

0.6176 

-0.1628 

0.31 

35 

0.9432 

0.6718 

0.5715 

-0.1410 

0.36 

(2) 

0.9433 

0.6717 

0.5714 

-0.1410 

0,37 

40 

0.7031 

0.5606 

0.5177 

-0.1211 

0.43 

(1) 

0,7032 

0,5607 

0.5177 

-0.1210 

-- 

(2) 

0.7029 

0,5605 

0,5177 

-0.1212 

0.43 

45 

0.5061 

0.4580 

0.4580 

-0.1026 

0.5 

(1) 

0,5063 

0.4582 

0.4582 

-0.1026 

0.5 

(2) 

0.5040 

0.4588 

0.4578 

-0.1024 

0 . 5 

(3)* 

0.5040 

0.4620 

0.4620 

-0.1026 

0 . 5 

*  (1)  Sattinger  and  Conway,  Reference  [10]. 

(2)  Morley,  Reference  [11]. 

(3)  Timoshenko  and  Woinowsky-Krieger,  Reference  [1]. 
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phenomenon  was  also  reported  by  Morley  [11].  The  boundary  errors  for 
the  deflection  and  slope  are  less  than  1.3%  and  2%,  respectively,  of 
the  maximum  values  of  the  deflection  and  slope  inside  the  plate  for 
a  >  30°.  As  the  angle  a  decreases  from  25°  to  10°,  the  boundary  errors 
around  the  acute  corner  gradually  increase.  Despite  this  difficulty  in 
the  vicinity  of  the  acute  corner  for  small  angles  a,  the  computed  num¬ 
erical  results  for  the  central  deflections,  bending  moments  and  maximum 
edge  moments  are  found  to  be  very  accurate  compared  with  the  results 

obtained  by  Morley  [11]  and  Conway  [10].  The  bending  moments  M  and 

/\ 

M  and  the  deflection  computed  at  the  plate  corner  are  plotted  in  Fig.  3. 

y 

The  numerical  results  after  satisfying  the  boundary  conditions 

r\  -  . 

w  =  0  and  -r—  at  8  and  10  equally  spaced  points  along  the  edge  AB  are 
listed  in  Table  3  together  with  the  results  obtained  by  9  points  matched 
along  edge  AB.  The  results  of  Table  3  indicate  the  convergence  of  the 
soluti on. 

3.3  Uniformly  Loaded  Circular  Segmental  Plates  with  All  Sides  Clamped 

Bending  of  the  uniformly  loaded  circular  segmental  plates  with 
all  edges  clamped  have  been  solved  by  Woinowsky-Krieger  [12]  using 
the  bipolar  coordinate  transformation  method.  An  alternate  approach 
to  solve  this  problem  using  point-matching  was  attempted.  The  sym¬ 
metry  of  the  problem  (see  Fig.  4)  shows  that  a  suitable  general  solu¬ 
tion  is 

w  =  D  [CT  (x4  +  2x2y2  +  y4)  +  C191  +  blGl 


■ 


•  it  gnafi 
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Table  3  Deflections  and  Bending  Moments  for  a 
Uniformly  Loaded  Rhombus  Plate  With 
All  Edges  Clamped,  Poisson's  Ratio  =  0,3 


Number 

a  of  points 

degrees  matched  on 
boundary 

Values  at  center  of  plate 

Max.  values  along  edge 

wxl02/(qa4/D) 

MxlO/qa2 

A 

MyXl0/qa2 

(M  )  /qa2 
v  n  max/M 

at  distance 
s / ( a/s i n  a) 

10° 

8 

3,0629 

1.3682 

0.70141 

-0,26005 

0.15 

9 

3,0213 

1.3596 

0.69212 

-0.26772 

0.14 

10 

2.9949 

1.3552 

0.68427 

-0.27265 

0.13 

15° 

8 

2.4495 

1.1981 

0.69819 

-0.24846 

0.18 

9 

2.4356 

1.1980 

0.69061 

-0.24992 

0.16 

10 

2,4278 

1.1982 

0.68609 

-0.24996 

0.15 

20° 

8 

1.9665 

1.0536 

0.67879 

-0.21999 

0.19 

9 

1.9638 

1.0543 

0.67662 

-0.21757 

0.18 

10 

1.9623 

1,0547 

0.67565 

-0.21455 

0.16 

25° 

8 

1.5673 

0.91887 

0.65266 

-0.18694 

0.22 

9 

1,5673 

0.91906 

0.65270 

-0.18465 

0.22 

10 

1.5676 

0.91915 

0,65276 

-0.18382 

0.23 

30° 

8 

1.2295 

0.79147 

0.61703 

-0.16175 

0.31 

9 

1.2298 

0.79145 

0.61728 

-0.16203 

0.31 

10 

1.2300 

0.79145 

0.61742 

-0.16217 

0,31 

35° 

8 

0.94316 

0.67191 

0.57144 

-0.14129 

0.36 

9 

0.94326 

0,67187 

0.57153 

-0,14097 

0.36 

10 

0.94331 

0.67185 

0.57156 

-0.14081 

0.37 

40° 

8 

0.70317 

0.56066 

0.51769 

-0,12111 

0.43 

9 

0.70316 

0 . 56065 

0.51769 

-0.12111 

0.43 

10 

0.70316 

0.56065 

0.51769 

-0.12118 

0.43 

45° 

8 

0.50613 

0,45809 

0.45809 

-0.10261 

0.50 

9 

0.506127 

0.45809 

0.45809 

-0.10261 

0.50 

10 

0.506125 

0.45809 

0.45809 

-0.10266 

0.50 

.  ' 

• 

y 

A 
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A 


(a)  a  >  j 


y 


(b)  a  <  y 


Fig.  4  Coordinate  system  for  circular  segmental  plate 


19 


N 

+k^l  ^ C4k-lg4k-l  +  c4kg4k  +  b4k-lG4k-l  +  b4kG4k^ 

Selecting  N  =  9  permits  satisfaction  of  the  boundary  conditions 
3  w 

w  =  0»  ^  =  0  at  19  points  along  AB  and  BC.  The  numerical  results  are 
shown  in  Table  4. 

The  locations  of  matching  points  along  the  boundaries  are  chosen 
as  follows: 

a)  In  the  cases  of  a  =  180°,  170°,  160°,  150°,  140°,  130°,  120°, 
110°,  100°  and  90°,  the  matching  points  are  so  located  that 
there  is  one  point  for  every  ten  degrees  along  arc  AB.  The 
rest  of  the  points  are  located  at  equally  spaced  intervals 
along  line  BC. 

b)  In  the  cases  of  a  =  80° ,  70°,  60°,  50°,  40°,  30° ,  20°  and  10°, 

3  w 

the  boundary  conditions  w  =  0,  =  0  are  satisfied  at  10 

equally  spaced  points  along  arc  AB.  These  boundary  conditions 
are  satisfied  at  9  equally  spaced  points  along  line  BC. 

The  numerical  results  are  shown  in  Table  4  and  compared  to  the 
results  from  Woi nowsky-Krieger  [12]. 

The  boundary  errors  for  the  deflection  and  slope  are  less  than 
0.1%  and  1%,  respectively,  compared  with  the  maximum  values  of  the  de¬ 
flection  and  slope  inside  the  plate  for  the  angles  a  >  20°. 

The  deflections  and  bending  moments  are  plotted  in  Fig.  5,  Fig.  6 
and  Fig.  7,  respectively.  The  maximum  moments  for  various  angles  a  are 
also  plotted  in  Fig.  8  where  the  dimensionless  quantity  is 
[M  /q(AC)2].  It  is  noted  that  the  reference  length  AC  is  used  rather 
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Table  4  Maximum  Deflections  and  Bending  Moments  for  a 
Uniformly  Loaded  Circular  Segmental  Plate  With 
All  Edges  Clamped,  Poisson's  Ratio  =  0.3. 


a 

degree 

Maximum 

deflect! ons 

(M  )  at  point  C,  M 

v  y'max  K  max 

,/qa 

wmax/(qa4/D) 

at  distance  s/AC 

This 

Work 

Woinowsky- 

Krieger 

10° 

0.1383 

X 

-9 

10  y 

0.50 

-0.1921 

X 

10-4 

-0.192 

x  10"4 

20° 

0.3409 

X 

10"7 

0.50 

-0.3018 

X 

10"  3 

-0.301 

x  10"3 

o 

O 

CO 

0.8191 

X 

10'6 

0.50 

-0.1481 

X 

10“  2 

-0.148 

x  10"2 

o 

o 

0.7469 

X 

10"5 

0.50 

-0.4475 

X 

10“  2 

-0.447 

x  10"2 

50° 

0.3949 

X 

10'4 

0.50 

-0.1029 

X 

10"1 

-0.103 

x  10"1 

60° 

0.1465 

X 

10"3 

0.50 

-0.1983 

X 

10"1 

-0.198 

x  10"1 

o 

O 

0.4217 

X 

_3 

10  J 

0.50 

-0.3361 

X 

10"1 

-0.337 

x  10“  1 

80° 

0.9893 

X 

10"3 

0.50 

-0.5161 

X 

10"1 

-0.519 

x  10"1 

90° 

0.2022 

X 

10"2 

0.52 

-0.7314 

X 

10"1 

-0.731 

x  10"1 

100° 

0.3582 

X 

10“2 

0.52 

-0.9675 

X 

10"1 

-0.964 

x  10'1 

110° 

0.5626 

X 

10'2 

0.52 

-0.1205 

-0.123 

120° 

0.8113 

X 

icf2 

0.52 

-0.1426 

-0.146 

130° 

0.1065 

X 

10'1 

0.52 

-0.1596 

-0.161 

140° 

0.1290 

X 

10'1 

0.52 

-0.1659 

-0.171 

150° 

0.1440 

X 

10_1 

0.52 

-0.1716 

-0.178 

160° 

0.1526 

X 

_  1 

10  1 

0.50 

-0.1767 

-0.170 

170° 

0.1559 

X 

icf1 

0.50 

-0.1582 

-0.154 

o 

O 

CO 

0.1562 

X 

10'1 

0.50 

-0.1250 

-0.125 
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Deflections  of  uniformly  loaded  segmental  clamped  plates 
along  the  center  line 
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Deflections  of  uniformly  loaded  segmental  clamped  plates  along 
the  center  line 
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CD 


Bending  moment  of  uniformly  loaded  segmental  clamped  plates,  along 
the  center  line 


24 


Fig 0  8  Maximum  bending  moments  for  uniformly  loaded 
segmental  clamped  plates 
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than  the  radius  a. 


3,4  Bending  of  Hydrostati cal ly  Loaded  Circular  Segmental  Plates 
with  All  Edges  Clamped 


Referring  to  Fig.  9,  all  the  conditions  for  this  problem  are  the 
same  as  in  the  above  section  (3.3)  with  the  exception  of  loading.  The 
hydrostatic  load  is  expressed  as; 

q  =  q0(l  -  £)  • 


The  particular  integral  to  be  used  is 


%  /X4  +  2x2y2  +  y4  x4y> 
T  ( - 64 - -  a24 


All  the  analyses  are  the  same  as  the  uniformly  loaded  circular  plate 
discussed  in  section  3.3.  The  final  numerical  results  after  satisfying 
the  boundary  conditions  are  shown  in  Table  5  for  a  =  40°  to  180°. 

When  the  angles  are  less  than  40°,  there  are  numerical  difficulties 
in  finding  the  coefficients  b.  and  c. .  Since  plates  with  small  angle 
a  are  not  usually  encountered  in  practice,  so  this  limitation  is  not  be¬ 
lieved  to  be  serious.  The  numerical  results  for  this  problem  is  be¬ 
lieved  to  be  new. 

The  boundary  errors  for  the  deflection  and  slope,  respectively, 
are  less  than  1.0%  and  3.4%  of  the  maximum  values  of  the  deflection  and 
slope  inside  the  plate  when  the  angles  are  a  >  40°. 
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Table  5  Maximum  Deflections  and  Bending  Moments  for  a 
Hydrostatically  Loaded  Circular  Segmental  Plate 
With  All  Edges  Clamped,  Poisson's  Ratio  =  0,3 


a 

degrees 

Maximum 

Deflect! on 

(Vmax.  at 

point  C 

a2) 

0 

w  /(a 
max7  VMo 

a  /D) 

at  distance  S/AC 

o 

O 

0 . 8822 

• 

10"6 

0.52 

-0,6318 

_3 

•  10  J 

50 

0  o 7083 

• 

10'5 

0.52 

-0.2221 

• 

o 

i 

)  PO 

60° 

0.3778 

• 

10“4 

0.52 

-0.6055 

•  10’* 

o 

O 

0o 1435 

• 

10"3 

0.52 

-0.1357 

•  10"1 

00 

o 

o 

0.4317 

• 

10“3 

0.52 

-0,2640 

•  10'1 

<£> 

o 

o 

0.1063 

• 

10~2 

0.55 

-0.4559 

•  10"1 

100° 

0.2213 

• 

10~2 

0.55 

-0.7126 

•  10"1 

110° 

0.3999 

• 

10"2 

0.55 

-0.1022 

o 

o 

C\J 

0.6431 

• 

10’* 

0.55 

-0.1365 

130° 

0.9267 

• 

10"2 

0.55 

-0.1688 

o 

O 

0.1202 

• 

10"1 

0.55 

-0.1901 

150° 

0.1406 

• 

10_1 

0.55 

-0.2097 

160° 

0.1531 

• 

10“ 1 

0.55 

-0.2267 

o 

O 

0.1578 

• 

10"1 

0.55 

-0.2088 

00 

o 

o 

0.1582 

• 

10 1 

0.55 

-0.1667 

Fig.  9 


Coordinate  system  for  hydrostatically  loaded  circular 
segmental  clamped  plate 


Fig.  10  Coordinate  system  for  uniformly  loaded  circular 
plates  with  diametrically  opposite  flat  sides 
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3.5  Bending  of  Uniformly  Loaded  Circular  Plates  with  Diametrically 
Opposite  Flat  Sides  with  All  Edges  Clamped 


Because  of  the  practical  importance  and  the  unavailability  of  the 
exact  solution  in  the  literature,  the  same  method  and  analysis  are  used 
to  obtain  numerical  results  for  this  problem.  Considering  the  symmetry 
of  the  problem  (see  Fig.  10),  equation  (5)  is  also  a  suitable  solution 
to  this  problem. 

The  numerical  results  after  satisfying  the  boundary  conditions 
w  =  °>  3n 

straight  edge  BC  at  a  total  of  10  equally  spaced  points  are  listed  in 
Table  6. 


9w  =  0  along  the  circular  arc  AB  and  w  =  0,  ~  =  0  along 


The  boundary  errors  for  the  deflection  and  slope,  respectively, 
are  less  than  0.1%  and  0.3%  of  the  maximum  values  of  the  deflection  and 
slope  inside  the  plate  when  the  angles  are  a  <  80°. 

The  numerical  results  of  clamped  rectangular  plates  with  the 
height,  2 ( A0 ) ,  and  width,  2(0C),  (see  Fig.  10)  are  also  given  in  Table 
6  for  comparison  with  the  results  of  uniformly  loaded  circular  plates 
with  diametrically  opposite  flat  sides. 


29 


Table  6  Maximum  Deflections  and  Bending  Moments  for  a 
Uniformly  Loaded  Clamped  Circular  Plate  With 
Diametrically  Opposite  Flat  Sides 


a 

w 

4 

/ qa  x  x=0 
'  D  '  y=0 

M 

ma^  at  C 
qa 

Mx 

2  x=0 

qa  y=o 

Note 

o 

O 

00 

0.3788  x  10~4 
*0.378  x  10'4 

-0.1005  x  10'1 
*-0.101  x  10"1 

0.5025  x  10'2 
*0.502  x  10‘2 

A0 

OC 

=  5.76 

o 

O 

0.5732  x  10‘3 
*0.573  x  10'3 

-0.3923  x  10"1 
*-0.392  x  10'1 

0.1962  x  10'1 
*0.196  x  10'1 

A0 

OC 

=  2.924 

O 

o 

uo 

0.2524  x  10"2 
*0.253  x  10'2 

-0.8271  x  10"1 
*-0.830  x  10'1 

0.4105  x  10'1 
*0.412  x  10"1 

A0 

OC 

=  2.0 

50° 

0.6045  x  10-2 
*0.614  x  10'2 

-0.1257 

*-0.127 

0.6121  x  10"1 
*0.617  x  10"1 

A0 

OC 

=  1.556 

o 

O 

0.1008  x  10"1 

-0.1547 

0.7433  x  10'1 

o 

O 

CO 

0.1316  x  10"1 

-0.1683 

0.7967  x  10"1 

20° 

0.1493  x  10'1 

-0.1637 

0.8113  x  10"1 

10° 

0.1556  x  10_1 

-0.1542 

0.8125  x  10'1 

0° 

0.15625  x  10'1 

-0.1250 

0.8125  x  10"1 

*  Results  from  clamped  rectangular  plate  with  height  2 (AO)  and  width  2(0C). 


CHAPTER  IV 


APPLICATIONS  TO  CREEPING  FLOW  PROBLEMS  IN 
SHALLOW  CAVITIES  WITH  VARIOUS  SHAPES 

4.1  Introduction 

The  motion  generated  in  a  fluid-filled  cavity  by  a  uniform  trans¬ 
lation  of  one  of  the  walls (see  Figs.  11,  12  and  18)  represents  one  of 
the  simplest  examples  of  the  steady  flow  involving  closed  streamlines, 
and  as  such  has  occupied  a  position  of  considerable  theoretical  import¬ 
ance  within  the  broader  field  of  steady  separated  flows.  Previous  work 
on  this  topic  has  been  reviewed  by  Burggraf  [13],  who,  for  the  special 
case  of  square  cavity,  also  obtained  numerical  solution  to  the  full 
Navier-Stokes  equation  for  a  range  of  Reynolds  number  R  =  VqD/v  from 
0  to  400  (V  is  the  velocity  of  the  top  wall,  D  is  the  width  of  the 
cavity  and  v  is  the  kinematic  viscosity  of  the  fluid). 

In  the  absence  of  the  inertia  terms,  the  equation  of  motion  reduces 

4 

to  the  familiar  bi harmonic  equation  V  ip  =  0,  where  ip  is  the  dimension¬ 
less  stream  function.  The  boundary  conditions  are: 

on  all  boundaries, 

at  the  upper  moving  wall, 

at  the  fixed  walls. 

Here,  n  is  the  outward  normal.  As  already  remarked  by  Burggraf  [13], 
the  form  of  the  boundary  conditions  precludes  an  analytic  solution  of 
this  system  by  one  of  the  standard  procedures  used  successfully  in  the 
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Fig.  11  Coordinate  system  for  a  rectangular  cavity 
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Fig.  12  Coordinate  system  for  a  triangular  cavity 
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field  of  elasticity,  or  by  approximate  methods,  such  as  variational 
technique. 

In  this  thesis,  the  point-matching  technique  was  applied  to  obtain 
approximate  numerical  results  for  cavities  of  rectangular,  triangular 
and  circular  segmental  cross  sections.  It  is  noted  that  the  particular 
integral  vanishes  for  the  creeping  flow  problem. 

4.2  Creeping  Flow  in  Cavities  of  Rectangular  Cross  Section 

Creeping  flow  solutions  for  rectangular  cavities  having  aspect 
H  11 

ratios,  A  =  —  =  1,  2,  3  and  5  were  obtained  numerically  by  a 

relaxation  technique  in  1966  by  Pan  and  Acrivos  [14].  In  this  thesis 
the  general  solution  of  biharmonic  equation,  g.  and  ,  in  Cartesian 
coordinates  was  used  to  obtain  an  approximate  numerical  solution  for 

l_l 

rectangular  cavities  with  aspect  ratio  A  =  —  =  1  and  2  by  the  point¬ 
matching  technique. 

The  symmetry  of  the  problem  (see  Fig.  11)  shows  that  equation  (6) 
is  a  suitable  general  solution.  As  noted  earlier,  the  particular  inte¬ 
gral  vanishes  for  the  creeping  flow  problem.  In  the  case  of  aspect  ratio, 

u 

A  =  —  =  1,  the  coefficients  c.  and  b.  are  determined  by  satisfying  ip  =  0 
a  i  i 

and  lit  =  1  at  5  equally  spaced  points  along  the  edge  CE,  ip  =  0  and 

lit  =  o  at  5  equally  spaced  points  along  the  edge  OB,  and  \p  =  0  and 
9y 

lit  =  o  at  9  equally  spaced  points  along  the  edge  BE.  The  numerical  re- 
suits  of  streamline  ip  and  velocity  distribution  along  the  center  line 
are  plotted  in  Fig.  13.  In  the  case  of  the  aspect  ratio  =  2,  the  numeri¬ 
cal  results  after  satisfying  the  boundary  conditions  \p  =  0  and  lit  =  1 
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Fig.  13  Creeping  flow  streamline  patterns  and  velocity 

distribution  along  the  center  line  for  square  cavity 


Cu  |  X 


34 


at  4  equally  spaced  points  along  the  edge  CE ,  ^  =  0  and  =  0  at  4 
equally  spaced  points  along  the  edge  OB,  and  ^  =  0  and  =  0  at  11 

o  X 

equally  spaced  points  along  the  edge  BE  are  plotted  in  Fig.  14. 

The  maximum  value  of  the  stream  function  at  the  first  vortex  is 
0.102  compared  with  the  value  0.1  given  by  both  references  [13,  14]  and 
the  location  of  the  vortex  center  is  -f-  =  0.235  compared  with  the  value 

a 

of  0.24  given  by  reference  [14]  and  the  value  of  0.25  given  by  reference 
[13]  for  a  square  cavity. 

The  maximum  and  minimum  values  of  the  stream  function  ip  for  the 
first  and  second  vortices  are  0.102  and  -0.000249,  respectively  (0.1 
and  -000023  by  reference  [14]).  The  vortex  center  locations  are 
|- =  0.26  and  1.55  respectively,  (0.25  and  1.575,  respectively,  by 

U 

reference  [14])  for  the  aspect  ratio  A  =  —  =  2. 

It  is  seen  that  the  results  from  the  present  analysis  are  very 
satisfactory.  The  only  disadvantage  of  the  point-matching  method  is 
that  the  lower  corner  vortices  which  have  been  discussed  by  Acrivos  [14] 
can  not  be  predicted  because  of  boundary  errors.  But  the  velocity  dis¬ 
tribution  along  the  center  line,  Y  axis,  can  be  obtained  easily  by  the 
point-matching  method.  This  is  a  very  important  property  of  the  problem. 
The  numerical  results  are  plotted  in  Fig.  13  and  Fig.  14  for  aspect 
ratios  one  and  two,  respectively. 

The  boundary  errors  for  stream  function  and  velocity  expressed  as 
a  percentage  of  their  maximum  values  are  less  than  0.6%  and  4%,  respec- 

u 

tively,  for  the  aspect  ratio  A  =  -  =  1.  The  boundary  errors  for 
stream  function  and  velocity  expressed  as  a  percentage  of  their  maximum 
values  are  less  than  1.0%  and  5%,  respectively,  for  the  aspect  ratio 


A  =  —  =  2 
A  a 


:  gw  f £ v  min  nr  i  t  5  it  I-xoh  !  ri 

* 


35 


Fig.  14  Creeping  flow  streamline  patterns  and  velocity  distribution 
along  the  center  line  for  rectangular  cavity  A  =  H/a  =  2. 
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4.3  Creeping  Flow  in  Cavities  of  Triangular  Cross  Section 


The  symmetry  of  the  problem  (see  Fig.  12)  shows  that  equation  (6) 
is  a  suitable  solution  for  the  triangular  cavities.  The  coefficients 
b..  and  c..  are  obtained  by  satisfying  ip  =  0  and  =  1  at  9  equally 
spaced  points  along  the  edge  AB  and  ^  =  0  and  =  0,  where  n  is  the 
outward  normal  along  the  boundary  BC,  at  10  equally  spaced  points  along 

i 

the  edge  BC.  The  numerical  results  of  the  stream  function  ip  and  the 
velocity  distribution  along  the  center  line  AC  are  plotted  in  Fig.  15, 
Fig.  16  and  Fig.  17  for  a  =  30°,  40°  and  50°,  respectively.  The  numeri¬ 
cal  results  of  streamline  \p  and  velocity  along  the  center  line  AC  are 
also  obtained  for  small  angles  a  =  10°  and  20°. 

The  numerical  results  for  iJj  M  and  its  location  are  listed  in 

max 

Table  7  for  angles  a  ■  10°  to  50°.  The  boundary  errors  for  stream 
function  and  velocity  expressed  as  a  percentage  of  their  maximum  values 
are  less  than  2 %  and  4.5%,  respectively,  for  angles  a  £  40°. 


4.4  Creeping  Flow  in  Cavities  of  Circular  Segmental  Cross  Section 

The  symmetry  of  the  problem  (see  Fig.  18)  shows  that  equation  (6) 
can  be  used  for  the  circular  segmental  cavity.  The  coefficients  b.. 
and  c.  are  determined  by  satisfying  ip  -  0  and  |y  =  1  at  9  equally 
spaced  points  along  the  edge  AB  and  ip  =  0  -  0  at  10  equally  spaced 

points  along  the  arc  BC.  The  numerical  results  of  the  stream  function 
\p  and  velocity  distribution  along  the  center  line  AC  are  plotted  in 
Fig.  19  and  Fig.  20  for  a  =  60°  and  90°,  respectively.  The  numerical 
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Fig.  17  Creeping  flow  streamline  patterns  and  velocity  distribution  along  the  center 
line  for  the  triangular  cavity  for  a  =  50° 
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Fig„  18  Coordinate  system  for  a  circular  segmental 
cavi ty 
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Fig.  19  Creeping  flow  streamline  patterns  and  velocity  distribution  along 
the  center  line  for  circular  segmental  cavity  for  a  =  60° 
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Creeping  flow  streamline  patterns  and  velocity  distribution  along 
the  center  line  for  the  circular  segmental  cavity  for  a  =  90° 
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Table  7  Numerically  determined  vortex 

center  location  in  triangular  cavity 


a 

ip 

ymax 

Vortex  Center 

(s/b) 

10° 

0.0126 

0.32 

20° 

0,0242 

0.30 

30° 

0.0350 

0.27 

o 

O 

■vt- 

0,0579 

0,23 

50° 

0,0663 

0,23 

Table  8 

Numerically  determined  vortex  center 
location  in  circular  segmental  cavity 

Vortex  Center 

a 

ip 

rmax 

(s/b) 

10° 

0.225 

X 

10"2 

0.35 

ro 

o 

o 

0.890 

X 

10"2 

0,32 

o 

O 

CO 

0.197 

X 

10'1 

0.32 

o 

o 

0.342 

X 

10_1 

0.32 

50° 

0.519 

X 

10'1 

0.32 

60° 

0.718 

X 

10"1 

0.32 

o 

O 

0.932 

X 

10"1 

0.32 

o 

o 

CO 

0.115 

0,32 

^ o 
o 

o 

0.136 

0,30 

100° 

0.155 

0.30 

110° 

0.170 

0.30 
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results  of  maximum  stream  function  ip  and  its  location  are  listed  in 
Table  8  for  angles  a  =  10°  to  110°. 

The  boundary  errors  for  stream  function  and  velocity  expressed  as 
a  percentage  of  their  maximum  values  are  less  than  1.5%  and  4%,  respec¬ 
tively,  for  a  <  60°.  As  the  angle  a  increases  from  80°  to  110°  the 
boundary  errors  around  the  upper  corner  gradually  increase.  The  maxi¬ 
mum  boundary  errors  for  velocity  is  about  7%  around  the  upper  corner 
for  a  =  70°. 

The  reason  that  the  large  boundary  errors  for  velocity  occur  in 
the  vicinity  of  the  upper  corner  for  all  the  creeping  flow  problems 
will  be  discussed  next. 

The  numerical  results  show  that  there  are  large  boundary  errors 

for  velocity  in  the  region  of  the  upper  corner  for  creeping  flow  problems 

due  to  the  discontinuity  of  the  velocity  in  this  corner.  It  is  physically 

impossible  that  a  velocity  field  is  discontinuous  at  a  certain  point 

since  the  shear  stress  should  be  infinite.  Mathematically,  the  solution 

of  bi harmonic  equation  should  be  continuous  in  the  second  derivative 

in  the  whole  region  of  the  problem.  However,  the  exact  solution  of  the 

creeping  flow  stream  function  subject  to  the  boundary  conditions  used 

must  have  a  discontinuity  in  its  first  derivatives  at  the  upper  corner 

E.  But  the  approximate  solution  is  written  in  terms  of  a  series  of 

q.  and  G.  which  are  continuous  in  their  first  derivative  (actually 
i 

continuous  up  to  nth  derivative,  where  n  can  be  any  positive  integer). 
Consequently,  the  error  in  the  first  derivative  of  the  approximate  solu¬ 
tion  is  large  in  the  neighbourhood  of  the  upper  corner. 
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CHAPTER  V 


CONCLUDING  REMARKS 


1.  The  exact  solution  in  Cartesian  coordinates  can  be  used  for  a 
class  of  problems  governed  by  the  bi harmonic  equation,  such  as 
uniformly  loaded  thin  plate  with  all  edges  clamped  and  creeping 
flow  problems. 

2.  The  well  known  general  solution  of  the  biharmonic  equations  in 
polar  coordinates  can  not  be  readily  applied  to  problems  such 

as  uniformly  loaded  clamped  rectangular  plates  with  large  aspect 
ratios  and  uniformly  loaded  clamped  rhombus  plates  with  small 
angles  a  (see  Fig.  1(a)  and  (b)).  In  contrast,  the  general  solu¬ 
tion  of  the  bi harmonic  equation  in  Cartesian  coordinates  yields 
the  approximate  solutions  to  the  above  problems  readily.  However, 
it  is  noted  that  the  expressions  for  the  harmonic  functions  g^ 
and  bi harmonic  functions  G^  become  quite  involved  as  the  number 
of  terms  to  be  used  increases. 

3.  The  numerical  results  indicate  that  the  point-matching  technique 

can  be  used  for  the  biharmonic  boundary  value  problems  when  the 

dW 

boundary  conditions  are  of  the  Cauchy  type,  i.e.  w  and  ^  are 
specified  along  the  boundary.  The  boundary  errors  are  very  small 
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compared  with  the  maximum  values  inside  the  region  for  most  of  the 
clamped  plates.  However,  for  the  clamped  rhombus  plate  with 
angles  a  less  than  25°  and  the  clamped  circular  segmental  plate 
with  angles  a  less  than  30°,  larger  boundary  deviations  were 
found. 

4.  Several  problems  with  boundary  conditions  other  than  the  Cauchy 
boundary  conditions  were  attempted  using  these  functions  g^  and 
and  the  point-matching  technique.  In  most  cases  boundary 
errors  were  large  in  comparison  to  the  maximum  values  inside  the 
region.  Reasonable  results  for  problems  with  non-Cauchy  type 
boundary  conditions  are  obtained  only  in  the  case  of  the  uniformly 
loaded,  simply  supported  rhombus  plate  with  angles  a  =  45°  and  40° 
(see  Fig.  1(b)).  Numerical  results  are  listed  in  Table  9,  where 
the  boundary  conditions  w  =  0  and  M  =  0,  are  satisfied  at  9 
equally  spaced  points  along  edge  AB.  The  boundary  errors  for 
deflection  and  bending  moment  are  less  than  0.4%  and  1.5%,  respec¬ 
tively,  compared  with  the  maximum  values  of  the  deflection  and 
bending  moment  inside  the  plate.  When  the  angle  a  is  35°,  the 
maximum  boundary  errors  for  the  deflection  and  bending  moment  are 
approximately  0.6%  and  7%,  respectively,  compared  to  the  maximum 
values  of  the  deflection  and  bending  moment  inside  the  plate. 
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It  is  also  found  that  errors  become  large  when  the  boundary  condition 
involves  the  second  derivative  of  the  deflection. 
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Table  9  Maximum  Deflection  and  Bending  Moments  of  a 
Uniformly  Loaded  Rhombus  Plate  With  All  Edges 
Simply  Supported 


degree 


45° 

0.16249  x 

i 

o 

o 

O 

0.22953  x 

10"1 

0.95771  x  10"1  0.95771  x  10 


0.11844 


0.10916 
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THE  G  FUNCTIONS 
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APPENDIX  2 
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=  0 


implying 


00  00 

I  I  a  (m) (m-1 ) (m-2) (m-3)  xm-4yn 
m=4  n=0  mn 


00  00 

+  2  l  l  a  (m)  (m-1 )  (n)  (n-1 )  xm"2  yn"2 
m=2  n=2  mn 


00  00 

*11  a  (n)(n-l)(n-2)(n-3)  xm  y11"4  =  0 

m=0  n=4  mn 


(Al) 


( A2) 


The  equation  (A2)  can  be  written  as: 


00  00 


implying 


[ak+4j,  (k+4) (k+3) (k+2) ( k+1 )  +  2ak+2j£+2  (k+2) ( k+1 ) 
U+2)U+1)  +  ak>£+4  (H+4)  (il+3)  (£+2) (£+1 )]  xky*  =  0 
[ak+4,il  (k+4)(k+3)(k+2)(k+l )  +  2ak+2j|l+2  (k+2) (k+1) 

{Jt+2)(«.+l )  +  ak  +4  U+4)(Jt+3)U+2)U+l)]  =  0 

whsi^G  k  _  0,l,2,3,.o..  00 
l  =  0,1 ,2 ,3,. ...  «> 


(A3) 


(A4) 
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Selecting  a  finite  number  of  terms,  then  equation  (Al)  can  be 


written  as: 


w, 


P  P 

z„  i 

m=0  n=0 


a  m  n 

am«  X  y 

mn  J 


(A5) 


where  m  +  n  <  p  and  p  >  4. 

Then  in  equation  (A4),  k  and  £  are 


k  =  0,  1,  2,  . 

I  =  0,  1,  2,  . 


but  k  +  £  <  (p-4) 


The  number  of  equations  (A4)  is  [ ( p- 3 ) (p-2)/2] 


For  instance,  p  =5 ,  w^  can  be  written  as: 


w 


a  +  a-,  +  a0>1y  +  a2>0x  +  a^xy  +  aQ>2y' 


h  a0,0  “1,0' 


3,ox  + 

a2,1x  y  +  *1,2^  +  *0.3* 

2  2 

3  . 

4 

2,2X  y 

+a1,3xy  +  a0,4*  +  a5,0 

3  2 

2  3 

4 

3,2X  ^ 

+  a2 ,3X  y  +  a4,ixy  +  a0 

5  4 

;  +  a4,lX  J 


(A6) 


Since  k  +  £  <  5  -  4  =  1,  then  there  are  three  constraint  equations 

(A4)  which  can  be  written  as: 

3a4 ,0  +  a2 ,2  +  a0 ,4  =  0 

5a5,0  +  a3,2  +  al ,4  =  0  ^ 

a4, 1  +  a2 ,3  +  5a0 ,5 


0 
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Omitting  all  the  details,  equation  (A6)  after  using  the  constraint 
equations  (A7)  can  be  written  as: 


w 


h 


I 

j=l 


11 


b.G.  + 
J  J 


i=l 


c.g 

iai 


where  the  coefficients  b.'s  and  c.'s  are  linear  combination  of  a  's 

j  i  mn 

and  all  the  coefficients  b.'s  and  c.'s  are  linearly  independent,, 
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